The relativistic glider revisited 
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In this paper we analyze some aspects of the "relativistic glider" proposed by Gueron and Mosna 
more in detail. In particular an explicit weak gravity and low velocity expansion is presented, 
the influence of different initial conditions are studied and the behavior of the glider over a longer 
integration time is presented. Our results confirm that the system can be used as a glider, but is 
not able to stop or even revert the fall of an object. 
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I. INTRODUCTION 

Changing the shape of a falling body in a gravitational 
field can induce changes in its motion. Such effects are 
well known and also could be used in space. In Newtonian 
gravity two body systems (two spacecraft connected by a 
tether) orbiting the Earth have been studied in [U 0, Q . 
Once-per-orbit modulations of the tether length [l|, yj or 
spinning tether systems with a modulation of the tether 
length linked to the spinning frequency [3( may be used 
as systems for propellentless propulsion. 

More recently Wisdom 0] and Gueron et al. stud- 
ied similar situations within general relativity. As most 
important difference the modification of the motion is no 
longer a resonant effect, thus the frequency of the change 
of shape is not linked to an orbital or spinning frequency. 
As pointed out in Q , their "gliding effect" has several ad- 
vantages over the "swimming effect" by Wisdom, in par- 
ticular it is not suppressed at high frequency and thus it 
appears to be possible to obtain a relevant displacement 
even in a very weak gravitational field such as the Earth 
more easily. However, the work of Gueron et al. @ only 
covers the radial motion and even there just presents the 
result of the integration over one period of cyclic changes. 
In this paper we present a more detailed analysis of the 
effect: it is identified in a weak gravity expansion, the de- 
pendence on the initial conditions is studied, it is shown 
how the system behaves over a longer integration time 
and finally we discuss different ways of implementing the 
system constraint. 



II. MODEL 

We consider the same model as introduced in Ref. [6|: 
two point masses connected by a massless strut which 
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moves in the gravitational field of Schwarzschild space- 
time. Analogously to that work, we use Schwarzschild 
coordinates as explicit coordinates and the motion is re- 
stricted to radial fall, such that the system is described by 
Schwarschild time t and the radii of the two point masses, 
rj. The massless strut is implemented as a time depen- 
dent constraint between the two radii: n(t) — ra(f)+i(t), 
where 



l(t;uj,a,Si) = l + Si exp 



(1 + a - 2utf 



(1 + a 2 )ujt{-l +uot) 



(1) 



for t € [0, T], and then is periodic. The time dependent 
length of the strut, l(t), is described by four parameters: 
its frequency w = 1/T, its amplitude Si, its minimum 
length and an asymetry parameter a. The asymmetry 
parameter a, taking values in the range [—1, 1], indicates 
how much the constraint fails to be symmetric with a = 
being the symmetric case. In Section I VII some issues 
about the implementation of the constraint are discussed. 

Putting the pieces together, the equations of motion 
are derived by using the action [Tlj 



S 



dt 



L 1 + ^L 2 + X(r 2 -r 1 -l(t)) 



(2) 



where A is a Lagrange multiplier enforcing the constraint 
and 



L i = l-Il-(l-Ii 



(3) 



with r s being the Schwarzschild radius. We suppose that 
the two masses are equal. In the following we also use the 
geometrical center of the system, r — (r\ + r2)/2. Then 
the system can be described with the variable r and the 
constraint l(t). 

We are interested in the difference between the vi- 
brating system and a non-vibrating system having the 
same initial conditions. This radial difference is noted 
5r(t) = r(t) — ro(t) where r(t) is the geometrical cen- 
tre of the vibrating system and r {t) is the geometrical 
centre of the non- vibrating one (our reference motion) . 
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FIG. 1: Comparison of the result from Ref. [g with the ex- 
pansion presented here, for R = 120 and r(t — 0) = 0. 



III. EXPANSION 

The initial radius of the apparatus is r(t = 0) = R. 
We first suppose that it is launched with a nul velocity 
r(t = 0) = 0. We introduce the order of magnitude 
O x = r s /R < 1. The quantity (Oi) Tl is denoted O n . 
We suppose that the amplitude of the constraint is small 
compare to the Schwarzschild radius, and is of the first 
order: Si/r s ~ 0\. Then we can expand the equations of 
motion, derived from the action @, around a reference 
motion: r = ro + Sr, where tq is the trajectory of the 
non- vibrating system, and Sr < Si. We suppose that we 
are in a weak gravitational field, and in the low velocity 
limit. Then we also expand the equations of motion with 
respect to r s /r and r. The resulting equations of motion 
are: 



Sr = 



31 
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The comparison between the exact result given in [6[ 
and the expansion is shown on Fig. [TJ The agreement 
is good for high oscillation frequencies. Indeed, the in- 
tegration time for each point is T = 1/lo. For low fre- 
quencies, the expansion is not valid anymore for two rea- 
sons: the distance covered by the system becomes non- 
negligeable compared to R, such that the r s /r expansion 
breaks down; the velocity of the system becomes non- 
negligeable compared to c, such that the low velocity 
limit breaks down. Additionally the validity of the ex- 
pansion imposes some limitations on the parameters u> 
and a. 

After one period T of oscillation, the velocity r'o ~ 
gT ~ Oi/uiR, where g — r s /R 2 , if the distance covered 
by the system during the time T is very little compare to 
R: r'oT -C R. Therefrom we obtain the condition ojR 3> 
0i/2, and we get tq <C Oi/2- On the other hand, Sr ~ 
LuSr < ujSi ~ bjRC>2- Supposing that loR < 0-3/2 > w e 
get Sr < 0i/2- The r s /r expansion is valid for Aro <C R, 
where Arg is the distance covered by the system during 




R (GM/c ) 

FIG. 2: Representation of the radial deviation 5r (obtained 
by solving the expanded equation of motion (Q) after one 
period of oscillation for ui — 0.05 and Si = 5 10 -3 . In the 
figure below, a log-log scale is used and in the top figure, a 
normal scale is used. The net displacment Sr is proportional 
to R~ 2 . This is confirmed here up to R = 10 11 . 



the time T. We have Ar - gT 2 - Oi/uj 2 R, thus leading 
to ujR 3> O1/2, which is coherent with the low velocity 
approximation. Then the expansion ([4]) is valid for 



(5) 



In Fig. [T]i? = 120 (GM/c 2 ) has been chosen, a value 
valid for a compact object such as a black hole or a neu- 
tron star. Then Oi ~ 1.7 ■ 10~ 2 and the expansion is 
valid for 1.1 ■ KT 3 < ujr s < 3.9 (in units of c 3 /GM). We 
can see from the condition ([5]) that for higher radius the 
expansion is valid for a broader range of frequencies. For 
Earth orbit with Oi ~ 10~ 9 , the expansion is valid for 
1.6 • 1CT 14 < u < 1.6 • 10 4 (in units of c 3 /GM). 

There is also a condition on a: if the assymetry pa- 
rameter a — > 1, then I — > 1 and the low velocity limit 
is no longer valid. Then the expansion is not valid 
when a is too near from one. 

In Q the net displacement after one period Sr = 
5r(T — l/u>) is estimated to be 



Sr ~ r(a 



Si GM 
R 2 ^ 



(6) 



where T is a dimensionless parameter that depends on 
the assymetry parameter a. In that work this relation 
is confirmed numerically for R ranging from 50 to 10 7 
(in units of GM/c 2 ). Using the expansion derived in this 
section we can confirm this relation for higher radii up 
to R = 10 11 . This in particular includes orbits around 
the Earth, where R ~ 10 9 . In Figure [5] the variation 
of Sr with repect to R as derived from the expanded 
equations Q are displayed. 



3 



IV. CHANGING THE INITIAL CONDITIONS 

In Ref. [6| the dependence of the displacement effect 
on the starting radius R was investigated and a l/R 2 
was found. However, no other dependence on the initial 
conditions was considered. 

Let us first look at different radial initial velocities, 
fit = 0) 7^ 0. In Figure [31 Sr for different initial veloc- 
ities at a positive value of a are presented. As can be 
seen, the effect increases if the initial velocity is negative 
(the system already falls down at t = 0), while in the 
opposite case (the system is thrown up at t = 0), the 
effect is decreased and at higher velocity becomes nega- 
tive. This behavior certainly is not surprising, since the 
work performed by the oscillations is expected to change 
sign when the system moves upwards. 
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FIG. 4: Representation of Sr as a function of u> (for R — 120 
and for a — —0.15) for different values of the initial radial 
velocity (noted dro expressed in fraction of c). For a fastly 
falling system it is still possible to get a positive effect. 
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FIG. 3: Representation of Sr as a function of u (and for 
a = 0.15) for different values of the initial radial velocity 
(noted dro expressed in fraction of c). We can see that for a 
falling system, the effect is increase while for a system going 
to upwards, the effect is decreased. 



a circular orbit around the central body which means 



GM 



(7) 



It can be seen that the effect decreases if the angular 
velocity is increased and becomes negative at high values. 
Still, the displacement is not very sensible to the initial 
conditions as long as the angular velocity remains small. 
This shows that the vibrations indeed can be used for 
"gliding," but will not help to increase the energy of an 
orbit. For further discussions on the effects of vibrating 
systems in an orbit around a central mass see Ref. 0]. 
Similar conclusions apply for negative values of a, but 
since for these values 5r < at — we do not reproduce 
these cases explicitly. 



Still, changing the sign of the asymmetry parameter 
for f(t = 0) > is not an option to obtain a positive dis- 
placement. This situation is depicted in Figure 2J where 
different initial velocities with negative values of a are 
presented. The change in the sign of a induces a vertical 
shift of all curves. Thus, a positive effect only is achiev- 
able for a system falling down with high velocity, while 
for initial conditions f = a negative value results over 
the whole frequency range, as was already found in Q. 

A different question concerns a non-trivial angular ve- 
locity as initial condition. For simplicity we assume that 
the strut stays aligned in radial direction during the com- 
plete motion. Thus it is sufficient to add one degree of 
freedom, the angle ip, to the system in order to allow ro- 
tations around the central body. This model is described 
more in detail in Ref. 0] • Figure [5] shows the behavior 
of Sr for different initial angular velocities. The angular 
velocities are represented in fractions of the angular ve- 
locity of a circular orbit. At (3 = the system falls down 
radially, at (3 = 1 the (non-vibrating) system is put into 
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FIG. 5: Sr as a function of lu for different values of initial 
angular velocity, expressed as a fraction of the circular angular 
velocity. In this example R = 120 and a = 0.4 has been 
chosen. 
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V. INTEGRATING OVER MANY PERIODS 

In Ref. [6| results have been presented from the inte- 
gration of one period of oscillation. This means on the 
one hand that the integration time is different for dif- 
ferent oscillation frequencies, on the other hand it is not 
obvious that the effect from the first period can be ex- 
trapolated to a longer integration time. As an important 
result Ref. @ obtained a characteristic "plateau" in the 
frequency dependence, as can seen in Fig. [TJ But since 
the integration time is not the same for different frequen- 
cies, one expect naively that this "plateau" represents in 
fact a linear increase of the effect with the frequency for 
a fixed integration time. 

Thus we present in Figure[5](5r as a function of time, in- 
tegrated over many oscillations. It can be seen that the 
displacement accumulates with time and for relatively 
small times, the expected linear increase for a fixed fre- 
quency is found. However, for low frequencies (ui = 0.02) 
the positive value of Sr turns into a negative displecement 
after some time, which is an effect of strong gravity field. 
Indeed, more simulations for different R have shown that 
the effect can increase with the time even for small fre- 
quency if R is sufficiently large. 
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FIG. 6: Representation of Sr as a function of the time (for 
R = 120, a = 0.15 and 5 t = 5 10" 3 ) for different values of u). 

These results show that for large initial radius R the 
naive extrapolation of the result of Ref. S] is correct. 
From the mentioned "plateau" (cf. Fig. []} an integration 
over a fixed time means that the total displacement grows 
linearly with the frequency. Thus, it is always interesting 
to push the frequency in order to enhance the total effect. 



VI. FERMI NORMAL COORDINATES 

As already mentioned in Ref. 6], the constraint is im- 
plemented in Schwarzschild coordinates and, physically, 
the strut should be constantly and locally monitored. 
This could be problematic as it appears to be more obvi- 
ous to implement the constraint by fabrication or selec- 
tion of a certain system, whose oscillation characteristics 



are measured in a laboratory beforehand. During the ra- 
dial fall experiment an active control of the oscillation 
should be avoided. In this situation the implementation 
of the constraint as a periodic motion in Schwarschild 
coordinates is incorrect. Instead, the constraint should 
have its prescribed shape in a coordinate system mov- 
ing together with the vibrating system. While possible 
in principle, the complexity of this calculation would be 
beyond the scope of this paper. On the other hand, 
we could implement the constraint in a locally inertial 
frame associated to the vibrating system. This is not the 
simplest method because one need to know a priori the 
trajectory of the vibrating system to describe its associ- 
ated locally inertial frame. We can also implement the 
constraint with a good approximation in the locally in- 
ertial frame associated to the non-vibrating system. Its 
trajectory is a free-fall motion (geodesic), which is non- 
accelerated (in the sense of General Relativity) and non- 
rotating (in the sense of Newtonian mechanics). The 
Fermi normal frame 8] is a realization of such a locally 
inertial frame. 




Constraint implemented in Swarzschild coord. 

Constraint implemented in Fermi coord. 
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FIG. 7: Representation of the constraint expressed 
in the Schwarzschild coordinates when implemented in 
Schwarzschild coordinates (dash-dot curve) and when imple- 
mented in Fermi normal coordinates (continuous curve). The 
constraint parameters are Si = 5 10 -3 , a = 0.4, lo = 0.05. 

The general expressions of the coordinate transforma- 
tions between arbitary coordinates and Fermi normal 
coordinates can be found in [t| flu ]. Using these ex- 
pressions, we calculate the transformations between the 
Schwarzschild coordinates and the Fermi normal coordi- 
nates associated to the radial free-fall up to fourth or- 
der. We implement the constraint ([TJ) in Fermi normal 
coordinates and apply the coordinate transformations to 
obtain the corresponding constraint in Schwarzschild co- 
ordinates. The comparison between the two different im- 
plementations of the constraint ([TJ) is shown in Figure [7] 
As expected, no essential difference can be seen after 
one oscillation period. However, the constraint imple- 
mented in Fermi normal coordinates fails to be periodic 
in Schwarschild coordinates and thus after many oscilla- 
tion periods the frequency between the two constraints 
is changed. Still, in a weak gravitation and over rather 



5 



short distances the influence of this shift is negligible. 

VII. CONCLUSIONS 

In this paper, we made a more detailed analysis of the 
relativistic glider, originally proposed by Gueron et al. 
■ We presented an analytical expansion that allows us 
to confirm the behaviour of the effect up to Earth 
radius. As a consequence, one can expect that it is not 
possible to see this effect in a very weak gravitational 
field. On the other hand, we showed that the devia- 
tion increases linearly with increasing integration time, 
which makes a long trajectory interesting for experimen- 
tal tests. 

We performed a more detailed analysis of the depen- 
dance of this effect on the initial conditions of the two- 
body system. If the system is already falling down when 
it starts to vibrate, the effect will be increased compared 
to the result of Ref. [6;]. On the other hand, if the sys- 
tem is thrown away from the central body, the effect 
decreases and eventually becomes negative. We also saw 



that a small tangential velocity has a very small impact 
on the effect, but it changes completely once the angular 
velocity is close to the value needed to get a circular orbit 
around the central body. For a more detailed study of 
this case we refer to Ref. [7|. 

Finally, it has been analyzed how the gliding effect 
changes if the constraint is implemented in the Fermi 
normal coordinates associated to the reference system in- 
stead of Schwarzschild coordinates. As expected, the re- 
sult does not change considerably in a weak gravitational 
field. 
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